Schwarzschild black hole as a grand canonical ensemble 



Gilad Gour* 

Racah Institute of Physics, Hebrew University of Jerusalem, 
Givat Ram, Jerusalem 91904, ISRAEL. 



Abstract 



For long black holes have been considered as endowed with a definite tem- 
perature. Yet when the Schwarzschild black hole is treated as a canonical 
ensemble three problems arise: incompatibility with the Hawking radiation, 
divergence of the partition function, and a formally negative mean-square 
fluctuation of the energy. We solve all three problems by considering the 
Schwarzschild black hole as a grand canonical ensemble, with the Hamilto- 
nian (the ADM mass) and the horizon surface area, separately, as observable 
parameters. The horizon area simulates the number of particles in statisti- 
cal mechanics since its spectrum is here assumed to be discrete and equally 
spaced. We obtain a logarithmic correction to the Bekenstein-Hawking en- 
tropy and a Gaussian type distribution for the energy levels. 



Typeset using REVTeX 



*E-mail: gour@cc.huji.ac.il 



1 



Black holes (BHs) are conventionally regarded as thermodynamic systems |T]||. But there 
exist problems with the description of a black hole as a canonical ensemble ||. For example, 
because a BH has negative specific heat, energy fluctuations calculated in canonical ensemble 
have formally negative variance. The issue of how to give a consistent thermodynamic 
description of a BH has received renewed urgency with the understanding that the mass 
spectrum of a BH may be discrete and highly degenerate. 

In the early seventies, Bekenstein pointed out that the horizon area of a non-extermal 
BH behaves as a classical adiabatic invariant Hf§]. Using Ehrenfest's principle |§ that 
any classical adiabatic invariant corresponds to a quantum entity with discrete spectrum, 
Bekenstein conjectured that the spectrum of the horizon area of a (non-extremal) BH should 
be quantized. He proposed a uniformly spaced spectrum which has later been considered by 
many authors (see the list in Kastrup @). Today the idea of a discrete eigenvalue spectrum 
for the horizon area is also supported by the work of Ashtekhar and others || , and recently a 
uniformly spaced spectrum (for large quantum numbers) has been established by Bojowald 
and Kastrup |5] within the framework of loop quantum gravity. 

Recently, the Schwarzschild BH with the uniformly spaced area spectrum was treated 
as a microcanonical ensemble by Scharf [JTTJ (a microcanonical ensemble of a gas of neutral 
and charged BHs has been considered earlier by Harms and Leblanc Jill) an d as a quantum 
canonical ensemble by Kastrup and by Makela and Repo |12|]. According to the "no 
hair" theorems |5|,|T3l|, an observer outside a Schwarzschild black hole is able to measure 
only its mass. Hence, all authors cited above have used the mass (energy) of the BH as the 
sole variable characterizing the system. The main new idea of the present paper is that a 
Schwarzschild BH formed by gravitational collapse should rather be considered as a grand 
canonical ensemble. The additional thermodynamic parameter is the horizon surface area 
of the BH. 

In the present paper we shall assume that the area spectrum A n of the Schwarzschild 
BH is given by: 
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A n = a n n = 0,1,2,... (1) 
where a is a constant of the order of the Planck area. Recently, Kastrup and others WA 



have shown that Eq. (|1|) holds true also when the space-time dimensions is taken to be 
D > 4. 

The area levels A n are expected to be degenerate. Denoting the degeneracy by g{n) and 
identifying lng(n) with the BH entropy 

S = A/(4C 2 p ) (2) 

(Cp and M.p will here denote the Planck length and mass, respectively), Bekenstein and 
Mukhanov |15| found that = 4£plnk, or equivalently g(n) = k n where k = 2,3,4,... 
They adopted k = 2 for simplicity; recently Hod [16j has argued in favor of the choice k = 3. 



In the quantum mechanical picture, if an observer at infinity makes an attempt to de- 
termine the mass of a black hole, his accuracy is limited by the time-energy uncertainty 
relation as well as by the systematic decrease of the mass of the BH |I7). Thus for a black 
hole formed by a gravitational collapse, it is impossible, even conceptually, for an observer 
to know exactly the mass of the black hole. However, if the black hole is in a static state 
(eternal black hole), the observer, in principle, would be able to determine its mass and then 
the microcanonical approach would be appropriate. 

In the canonical ensemble approach the observer knows only the average value of the 
mass of the black hole at a given time. In this ensemble the partition function Z formally 
diverges because the degeneracy factor g{n) rises very fast with n. Kastrup ]7j proposes to 
resolve this problem by defining Z by analytic continuation, a procedure which is difficult 
to understand in physical terms. Makela and Repo ]12| avoided the problem by studying, 



not the partition function of the whole spacetime itself, but instead the partition function 
of the radiation emitted by the BH. We shall see shortly that if we treat the BH as a grand 
canonical ensemble, Z is no longer divergent! Before describing this new approach we shall 
first point out another problem of the canonical approach. 
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Treating the BH as a canonical ensemble implies that it is in a thermal equilibrium with a 
surrounding thermal bath. Thus the BH mass remains constant because the density matrix 
of the canonical ensemble is constant in time; it commutes with the Hamiltonian (the ADM 
mass operator). However, the mass of an isolated BH formed by a gravitational collapse 
decreases in time because of losses to Hawking radiation [F7|], so a description of it via a 



thermal ensemble seems inappropriate. We propose to solve both problems by abandoning 
the canonical approach in favor of the grand canonical one. 

If the observer is interested in determining the horizon area of the black hole he is limited, 
apart from the time-energy uncertainty relation, by some kind of area-phase uncertainty 



relation 18 



AAA0 > ~a Q . (3) 

Here is the canonical conjugate to the area (number) operator. This is the first clue 
suggesting to treat the black hole as a grand canonical ensemble: the observer only knows 
the average value of the horizon area (by Eq. (|3|) the area operator is a number operator, and 
is thus analogous to the number of particles in a grand canonical system). We shall discuss 
now further physical grounds for adopting the grand canonical approach for the description 
of Schwarzschild BHs. 

Classically, the ADM mass and the horizon surface area are related by 

A = - — M . (4) 

c z 

Hence, the horizon area was never before considered as a new parameter. However, the mass 
(energy) and the horizon area (topology feature) of a BH describe two different properties. 

Consider an observer who is not aware of the classical relation in Eq. (|4]), and is interested 
in the BH properties. It is clear that he would use completely different techniques (and 
apparatus) to measure the mass and the horizon area of the BH. The mass may be measured 
asymptotically at infinity whereas the horizon surface area should be measured locally. Thus, 
these two parameters describe operationally distinguishable features of the BH. 



In the quantum mechanical picture the distinction between these two parameters becomes 
prominent. Because of the Hawking phenomenon, the horizon area decreases in time and 
hence the area operator and the Hamiltonian (the ADM mass operator) do not commute. 
Thus, the relation in Eq. does not hold true for operators. In |18j we have found the form 
of the Hamiltonian and have shown that in the classical limit, where the HR is negligible, 
Eq. (U) is satisfied also for the operators. 

Generally, in grand canonical ensembles the number of particles and the energy of the 
system are taken to be the two observable parameters since neither is constant and each 
describes a different property of the system (even though they are ultimately related by some 
formula). Since the ADM mass and the horizon area of a Schwarzschild BH have the same 
relation, we conclude that a Schwarzschild BH should be considered as a grand canonical 
ensemble! 

Assuming the black hole is described by some density operator p, we shall maximize the 
following quantity (entropy): 

Q = -Tr (p lnp) — p'(A) — b(H) (5) 

where p' and b are the Lagrangian multipliers (the physical meaning of these parameters 
will be discussed later) and H is the Hamiltonian (ADM mass) operator (boldface is used 
everywhere to denote operators). The trace in Eq. @ may be taken with respect to the 
area eigenstates; this makes it easy to take the degeneracy factor g{n) into account. The 
extremum for Q under the conditions that Tr (pA) = (A) and Tr (pH) = (H) with (A) and 
(H) known is attained by 

p = iexp(-p'A-6H), (6) 

where the partition function Z is defined by 

Z = Tr (exp(-p'A - 6H)) . (7) 

Because of the Hawking radiation, the Hamiltonian operator cannot commute with the 
horizon area operator, because the last operator is not constant in time. Thus p does 
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not commute with H. Hence, we conclude that by choosing the appropriate Hamiltonian 



operator [pl| for the Schwarzschild BH, the density matrix in Eq. (10) would be compatible 



with BH evolution in the wake of Hawking radiation. 



We have shown in ref. ]Tj| that the Hamiltonian (the ADM mass) can be written as 

H = M + V (8) 



where M = yc 4 A/167rG 2 is the mass operator (according to the classical limit) and V is a 
coupling between the horizon area and its canonical conjugate, the phase of the BH. However, 
as we have pointed out, the interaction term approaches zero as (M) — > oo according to 
(V) ~ A4pC 2 / (M). Hence, for BHs not near the Planck scales we may neglect the effect of 
Hawking radiation and assume that H « M. Thus, the partition function given in Eq. (^) 
can be approximated by 

oo oo 

Z = ^2 k n exp (^—fi'a n — bm y/rij = exp (/i'a — In k)n — bm y/n^j , (9) 

n=0 n=0 



where m = J c 4 ao/16irG 2 is of the order of the Planck mass. Let us define the dimensionless 
coefficients a = aofi' — In k and x = bmo. We note that the partition function converges for 
a > 0; we thus assume a > and that a is crudely of order unity. The probability to find 
the system in the nth area state can be written as 

P n = \ exp (-an - xV™) (10) 



Z 

where a and x must be reexpressed in terms of (A) and (M). 

The relation between (A) and (M) is not exactly (A)=167r(M) 2 G 2 /c 4 since the expec- 
tation values are not taken with respect to a pure state with sharp mass. What would it 
be for a BH formed by gravitational collapse? In order to answer, let us first introduce a 
new parameter n , the number n that maximizes P n . It is clear that around this value the 
distribution is symmetric, that is 

Pno+h ~ Pno-h (H) 
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for /i < no. Comparing Eq. ( pTTj) with Eq. (|T0|) we find that 

X = -2y^a. (12) 

Eq. flT2| ) raises a serious problem regarding the physical meaning of the parameters // and 
b. A comparison of the partition function ([5]) with the partition function 

Z = J2 g(E n ) exp (J3(jiN n - £„)) (13) 

n 

written for a general grand canonical ensemble with energy levels E n , particle number N n , 
degeneracy g(E n ), chemical potential \i and inverse temperature (3 reveals that the param- 
eter y! in our analysis represents the negative of the chemical potential divided by the 
temperature, and parameter b the inverse temperature of the BH. Now, Eq. (|T2|) implies 
that x is negative because a > and thus also b = x/ m o is negative. Does this imply that 
either the BH's mass or its temperature are negative? 

In the view of Makela and Repo [f[^| ; one can solve the problem by assuming mo is 
negative. In that case, b is positive (and thus the temperature too) and M n increases when 
n decreases. In other words, M n becomes greater when the BH becomes smaller. The 
meaning of M n changes: it is not the mass of the BH but, the energy of the BH radiation 
(with an appropriate choice of energy zero) . Here we suggest another solution to the problem 
which also saves the positivity of both the mass and the temperature of the BH, even though 
b is negative. 

Since the horizon splits the space into two parts, we cannot immediately compare our 
analysis with the one in general grand canonical ensemble. According to the "no hair" 
theorems, the degrees of freedom in the interior region are not accessible to an observer 
at infinity. This affects the independence of the horizon surface area and the mass of the 
BH. Independence of A on M implies, for example, that it is possible to change A (at least 
slightly) while keeping M constant. Hence, at infinity one cannot observe the independence 
of A and M due to the "no hair" theorems. 

The partition function that maximizes the entropy of the whole BH spacetime is given 
by Eq. @. Now, in order to associate with the BH a temperature Tbh = h/8nM, we 
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must restrict ourselves to the exterior region of the BH; there is no meaning to temperature 
"inside the BH" . How is this restriction implemented in our scheme? 

According to the above arguments we have to make the transition from two independent 
parameters A and M to one parameter. Normally, in the usual grand canonical ensemble, 
the parameter b would be given by 



b =[m) w (14) 

Note that this derivative would be taken with respect to (M) while keeping (A) constant, 
and would be negative. But when we are restricted to the exterior region, there is only 
one parameter, say (M) , and therefore the inverse temperature is defined by (the entropy is 
derived in Eq. ( p0|) ) 

. f dS \ 8tt(M) / 1 \ 

which is positive and equal to the inverse of Hawking's temperature. Because the derivative 
in Eq. fllBD is the total derivative, it is distinct from b, and can be positive. No negative 
temperature is necessary. 

Substituting Eq. (|TJ) back in Eq. ([TOP we find in the limit n — > oo 

P n = 7?exp (-a(n - « (^-) ex P ( _ 4~(™ ~ ™°) 2 ) ■ ( 16 ) 



Hence, we obtain a Gaussian distribution with a variance a a = y 2ri /a. Note that n is 
approximately the average of N = A/ao- Thus, as is typical of many-particle statistical 
systems, as (N) — > oo, the absolute fluctuations become large, but the relative fluctuations 
approach zero. 

The entropy may now be expressed as 

S = -Tr (plnp) = ii'(A) +b(M) +\nZ (17) 

where for large no the partition function || is given up to a very good approximation by 



Z = exp(ano) W u - + 0(l/n^) . (18) 



Calculating (A), (M) with (|16|) and taking into account the first order corrections we find 

(A) = a n + — ^ + O(l/n ) 

(M) = m 0v ^ + + (l/nf ). (19) 



Substituting all these in Eq. (|T7|) we finally obtain 

1 ... 1 /(A)\ 1 /4tt\ 3(lnfc-l) 1 
g = 77^ A > + oMW +oM~ + - (20) 



4£2 N ' 2 V y 2 U/ 2a 2 ' 

Note that the main contribution to the entropy is given by (A)/4C 2 as was to be expected. 
The logarithmic correction to the entropy is exactly the same as Makela and Repo obtained 



for the emitted radiation [12] and as Kastrup obtained from his analytic continuation ap- 



proach 0. However, our grand canonical approach does not suffer from a divergent partition 



function. Furthermore, using Eq. (19) to determine the fluctuation in the mass we find 



al = (M 2 ) - (M) 2 = (21) 

ana 

Thus, the mean square fluctuations of the energy is positive even while the specific heat is 
negative. 

In summary the grand canonical approach solves three problems which arise when us- 
ing the canonical approach with the area spectrum given in Eq. ([!]). Firstly, the partition 
function is automatically convergent. Secondly, the grand canonical approach is compatible 
with Hawking radiation in the sense that it requires the density operator to vary with time. 
Thirdly, the mean square fluctuations of the energy comes out positive. The distribution 
of the area (energy) levels is of the Gaussian type (for (A) £ 2 ), with relative fluctua- 



tions AA/(A) of order of C P /J(A). Furthermore, as a byproduct of the grand canonical 
approach we have recovered the same logarithmic correction to the entropy earlier derived 



by Kastrup |7j and by Makela and Repo ||12|| . Other authors |19|-t2T|| have also obtained a 
logarithmic corrections to the entropy of BHs. 
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